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Abstract 

Based on the two-dimensional lattice fermion model, we discuss transitions 
between different pairing states. Each phase is labeled by an integer which is a 
topological invariant and characterized by vortices of the Bloch wavefunction. 
The transitions between phases with different integers obey a selection rule. 
Basic properties of the edge states are revealed. They reflect the topological 
character of the bulk. Transitions driven by randomness are also discussed 
numerically. 
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Quantum phase transitions between different superconducting states have attracted much 
interest recently. In refs. [|l],0, for example, its possible realization in a high-Tc supercon- 
ductor was proposed, which is accompanied by the time-reversal symmetry T breaking. 
Further, there is a recent observation that it has some similarity to the plateaux transition 
in the integer quantum Hall effect (IQHE) One of the claims is that each phase is 

labeled by an integer (an analogue of the Hall conductance in the IQHE) and there can be 
transitions between phases with different integers. 

In this paper, based on the lattice fermion model, we investigate the problem. The 
integer for each phase is defined by a topological invariant of the U(l) fiber bundle (the 
Chern number) The U(l) fiber bundle is a geometrical object which is composed 

of the Brillouin zone (torus) and the Bloch wavefunctions (fiber). Due to its topological 
stability, a singularity in the U(l) fiber bundle necessarily occurs with the change of the 
Chern number. The singularity is identified with the energy-gap closing [pHT^. The Chern 
number is closely related to the zero points (vortices) of the Bloch wavefunction. Focusing 
on the motion of the vortices near the singularity, we give a general proof of a selection rule 
of the transitions. Due to the intrinsic symmetry of the system, the selection rule differs from 
that of the IQHE . We also investigate the properties of the edge states and how they 
reflect the topological character of the bulk. The transition due to the change of randomness 



strength is a typical example of the problem. As emphasized in ref. |T^, the symmetry effect 
leads to a new universality class and it is interesting as the Anderson localization problem. 
We also discuss the disorder-driven transition numerically. 
The Hamiltonian is 

l,m l,m ^^ml ^rnl / 

where cj^ = (c]^|,cj^|), c„ = *(c„|,c„|) and n = (ra^., ny)GZ^. This is an extension of the 
lattice fermion model discussed in connection with the plateaux transition in the IQHE 
||TO|-p!^. Here we comment on the relation between this Hamiltonian and the superconduc- 
tivity. Under the unitary transformation c„|— c„|— (for Vn), the Hamiltonian ([l|) is 

equivalent to H = J2i,m[d^timdm^ + d\^timdra[ + d\^/\imd^^^ + dm[^*imdi'\\- This is the pairing 
model for the singlet superconductivity. In the context of the superconductivity, the pair 
potential A^^ should be determined by the self-consistent equation. Although the effect is 
interesting in itself, it is beyond the scope of this paper. Further, the conditions t;*^ = t^i 
and A;m = ^mi are imposed and they correspond to the hermiticity and the SU(2) symmetry 
respectively. The SU(2) symmetry leads to the condition 



- a. 



Himay)* = 'Rim- (2) 



Due to the SU(2) symmetry, we can restrict ourselves to the sector d]^o"^d„ = without 
loss of generality. This is equivalent to the half-filled condition for the Hamiltonian (P, 
which we impose in the following arguments. 

Now let us define a topological invariant (the Chern number) for our model. It is a key 
concept in the following arguments. Put the system on a torus, which is x Ly and periodic 

in both X and y directions. Define the Fourier transformation by c„ = 1/ ^JlI^X^J2k^^'^"^^{k) 
where k = [k^, ky) is on the Brillouin zone (— vr, 7r]x (— tt, tt]. Assuming to be invariant 
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under translations, we obtain H = X]a:C(A;)'''H(A;)c(A;) where H(A;) = I](;_m)6~*^^'~'"^H;m- 
The H(/c) has two eigenvectors and eigenvalues. They correspond to the Bloch wavefunctions 
and the energy bands respectively. To satisfy the half-filled condition, the lower band is 
occupied for the ground state. We denote the Bloch wavefunction for the lower band by 
* {a{k),b{k) ). Then the topological invariant (the Chern number of the U(l) fibre bundle) 
is defined as: 

C=^ fdkz-{VkxA) (3) 
Zm J 

where A = {a*{k),b*{k))Vk^{a{k),b{k)) and z = (0,0,1) The integration J dk is 

over the Brillouin zone which can be identified with a torus. For simplicity, we assumed 
Him to be invariant under translations. However a generalization to a multi-band system 
(including a random system), is possible. It is crucial for the following arguments to rewrite 
the above formula in terms of a zero point of the Bloch wavefunction {vortex) and the winding 
number (charge). To be explicit, let us perform the gauge fixing of the Bloch wavefunction 
for the lower energy band. We note that the Chern number itself does not depend on the 
gauge fixing. To define the gauge, we use the rule a{k) = 1 and introduce a notation 
b{k) = b'{k)e~^'^^''^ {b'{k)EH). An ambiguity in the gauge fixing occurs when a{k) = 0. 
Around the zero point (vortex) in the Brillouin zone, it is necessary to change the way of 
the gauge fixing, for example, as b{k) = 1. Then the Chern number (^ is rewritten as 

C = Y.Ci, Ci = 1/271 i dkVCik) (4) 

where the summation is taken over all vortices of a{k) and Ri is a region surrounding the 
l-th vortex which does not contain other zeros of either a{k) or b{k). Here Ci is an integer 
and we call it the charge of the l-th vortex. 

Let us discuss the dx2_y2 + idxy model on a torus as an example The model is 

defined by tn+ea;,™ tn+ey,n ^n+ex,n ^n+Cy.n ^x^—y^^ ^n+ea;+ej;,n ^n—ex-\-ey,n 

il^xy for (e^. = (1,0), = (0,1), t > 0, A^.2_j^2, A^,j^gR) and the other matrix elements 

/ M]A B(k) \ 

are zero. Then H(A;) = {^j^^^j^-^ A{k)) ^(^) — 2t(cosA;^. + cosA;^^) and B{k) = 

2Aa.2_y2(cos fcj: — cos ky) +2iAxy[cos {kx + ky) — cos {k^ — ky)]. The energy spectrum is given 



hj E = ±yA(A;)^ + \B{k)\'^. When A^.^ = 0, the upper band and the lower band touch 
at four points (±7r/2, ±7r/2) in the Brillouin zone. The low- lying excitations around the 
gap-closing points are described by massless Dirac fermions. By turning on a finite ^xy, 
the mass generation occurs in the Dirac fermions. The vortex position is given by a{k) = 
and it is (k^^ky) = (0,0). Using B{k) = B\k)e'^^^^ (5'(fc)GR), the charge of the vortex is 
l/27r/(Q Q-) dkVQ{k) = 2 sgR{Axy/ Ax2^y2). Since there is no other vortex, the Chern number 
C is given by 

C = 2 sgn{Axy/Ax2_y2). (5) 

Here we note that, as in the case of the IQHE on the lattice, the Chern number can take 
various integer values in general cases e.g. a multi-band system and a model with a different 
T-broken pairing symmetry. 
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As in the QHE [p!^^ -[T9|, the edge states play a crucial role in the problem 0-^. The edge 
states reflect the bulk properties and it is possible to detect the topological character of the 
bulk through the edge states. In order to discuss the edge states, put the system on a cylinder 
which is Lj-xLy and periodic only in y direction. Further we impose an open boundary 

condition in x direction. Define the Fourier transformation by c„ = 1 / ^Jl^J^ky^^''^^^ '^n^iky) 
where ky is on (— vr, tt]. Then H = J2kyCia:{kyyiii^m^{ky)Cm^{ky) where Him is assumed to be 
invariant under translations in y direction and iii^m^{ky) = I](/y-my)6~*'^^*-'''~™'^''H/m- The 
relation (0) can be rewritten as —{ay'H.i^rn^{ky)ay)* = lii^m^{—ky). Define an eigenvector 
u by J2m^ ^ixm^{K)'^m^ = ^Ui^. Now wc show that there are two basic operations V and 
Q on the vector. They are defined by ('Pu)„^ = {(JyUnJ* and (Qu)„^ = ul^.^^+i. Then 
H/^m^(-^j/)('Pu)m^ = {-E){Vu.)i^. Further, when tim is real and uniform, we can 
obtain another relation Y^mxiii^m^{—ky){Qu)m^ = E{Qu)i^. Based on the symmetry, we 
shall discuss basic properties of the edge states. Consider a case when u is an eigenstate 
which is localized spatially on the left (or right) boundary i.e. a left (or right)-hand edge 
state. Then, from the above argument, u,Vu,Qu,V Qu are classified into two left-hand edge 
states and two right-hand edge states. Now, as in the argument of the IQHE ||2^, let us 
introduce a fictitious flux through the cylinder and change it from to flux quanta hc/e. 
Due to the symmetry, the number of edge states which move from one boundary to the other 
is necessarily even. We shall consider the dx-2-y2 +idxy model on a cylinder as an example. In 
Fig.l, the energy spectrum is shown. It can be confirmed that the energy spectra of the edge 
states appear in pairs and the number of the edge states which move from one boundary to 
the other as the fictitious flux is added, is even and coincides with the Chern number. In 
other words, the edge states directly reflect the topological character of the bulk. 

As discussed above, each phase in our model is labeled by the Chern number. The vortices 
move as the Hamiltonian is perturbed. The Chern number, however, does not change in 
general. Due to the topological stability, the change of the Chern number is necessarily 
accompanied by a singularity in the U(l) fiber bundle. The singularity is identified with the 
energy-gap closing. Focusing on the singularity, we can prove a selection rule from a general 
point of view (see also [pUj-p!^). The selection rule is closely related to the SU(2) symmetry in 
our model. Let us introduce a parameter g in the Hamiltonian. Assume that, when g = go, 
the energy gap closes at several zero-energy points in the Brillouin zone. Next focus on the 
region near one of the gap-closing points (/c°, k^) i.e. p = ^ {k^ — k^, ky — ky, g — go) ~ 0. 
Then the leading part of the Hamiltonian is, generally, given by 

Ho(p) = Ivop + (cr^, ay, cr^)vp (6) 

where vq is a 1x3 vector, (Jx{y,z) is a 2x2 Pauli matrix and v is a 3x3 matrix. Now 
let us introduce the standard form, which is convenient for the following arguments 
rT| . Choosing a unitary transformation U appropriately, one can obtain WHo(p)W~^ = 



IvqP + {ax, ay, (T^)DTp where D is diag(l, 1, sgn(detv)) and T is an upper triangle matrix 
with positive diagonal elements. Let us perform Tp p (the parity-conserving Affine 
transformation on [kx, ky) and the rescaling on g) and the redefinition voT^^^vq. Finally 
the standard form is obtained as 

Hi(p) = Ivop + axPx + ayPy + cr^p2Sgn(det(v)). (7) 
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This is equivalent to the Hamiltonian ii{k) where A{k) = p2Sgn(detv) and B{k) = — ipy. 
Performing the same analysis, one can find that a vortex moves from one band to the 
other at the gap-closing, and the conclusion is that the change of the Chern number is 



practically determined by sgn(detv) and the change is +1 or —1 |T^. Next let us consider 
a dual gap-closing point {—k^, —ky) which exists due to the symmetry. Here the relation 
derived from plays a crucial role and it is given by H(— A;) = —{o'yli{k)o'y)* . Therefore 

H(-fc° - p^, -/c° - py, g) = -(cryH(fc° + p^, + py, q) a yf Ivop + (cTx-, (Ty, cr^)vp. To 

summarize, the linearized Hamiltonian near the dual gap-closing point {—k^^, ~^y) given 
by 

H(-A;° + p^, -kl + py, g)r^ - IvqP + (cr^., cr^, cr^)wp (8) 

where w = vdiag(— 1, —1, 1). It gives sgn(detv) = sgn(detw). Therefore the change of the 
Chern number due to gap-closings always occurs in pair with the same sign and the total 
change is AC = ±2 generally. This is the selection rule. On the other hand, in the absence 
of the relation (j^) (or SU(2) symmetry), the above argument does not hold and it leads 
to the rule AC = ±1. Now we note the results in ref. 0] where the network model with 
the same symmetry as our model was investigated. In spite of the fact that their model is 
different from the lattice fermion model considered here, our selection rule still applies: this 
suggests the universality. Assuming that the system can belong to a phase with a vanishing 
Chern number by tuning parameters in the Hamiltonian, the selection rule implies that the 
Chern number is always even, which supports the result based on the edge states. 

Finally, we comment on the disorder-driven transition based on the random d^2_y2 +idxy 
model. It is defined by tij = t^- + 5tij and Ajj = A°. + 5/S.ij where 5tij and 5Ajj denote 
the randomness. Here the hermiticity and the SU(2) symmetry are imposed on tij and Ajj 



respectively. It has an intimate connection with the random Dirac fermion problem |21-26 



It is to be noted that the SU(2) symmetry is kept even in the presence of randomness and 



the model is interesting as the Anderson localization problem [|T6|. As discussed above, the 
Chern number is ±2 in the absence of randomness. On the other hand, in the presence of 
sufficiently strong randomness, it is expected that all the vortices disappear through the pair 
annihilation of vortices with opposite charges and the Chern number vanishes 0,0]. By 
the numerical diagonalization, we treated the disorder- driven transition for the 5tij = 6ijfi, 
SAij = SijQi where the f^s,g^s are uniform random numbers chosen from [—W/2, W/2]. The 
model was also studied extensively in ref. [^]. In Fig.2, the density of states is shown in 
the case A^^t^O. It can be seen that the two energy bands become closer and finally touch, 
as the randomness strength is increased. The transition C = ±2^0 with the gap-closing is 
a natural consequence from the selection rule. The exploration of the global phase diagram 
and the field-theoretical description are left as future problems. 

The authors thank T. Fukui for valuable discussions and P.-A. Bares for careful reading 
of the manuscript. This work was supported in part by Grant-in- Aid from the Ministry of 
Education, Science and Culture of Japan and also Kawakami Memorial Foundation. The 
computation has been partly done using the facilities of the Supercomputer Center, ISSP, 
University of Tokyo. 
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FIGURES 



FIG. 1. Energy spectrum in the d^2_y2 + idxy model on a cylinder {t = 1, Aj.2_j^2 = 1 and 
= 0.5). The system size is 40x40. The spectrum in the energy gap corresponds to the edge 
states. The boundary where the edge states are localized is shown by L(left) or R(right). 

FIG. 2. Density of states in the random dr^2_y2 + idxy model. In the absence of randomness, 
the model reduces to the pure dx2_y2 + idxy model with t = 1, Ax2_y2 = 1, and Axy = 0.5. 
The system size is 30x30 and the ensemble average is performed over 800 different realizations of 
randomness. 
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